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Abstract. Multi-dimensional constrained covariant density functional theories were developed recently. In 
these theories, all shape degrees of freedom 6,, deformations with even u are allowed, e.g., B29, B22, B30, B32, 
Bao, B42, B44, and so on and the CDFT functional can be one of the following four forms: the meson exchange or 
point-coupling nucleon interactions combined with the non-linear or density-dependent couplings. In this contri- 
bution, some applications of these theories are presented. The potential energy surfaces of actinide nuclei in the 
(B20, B22, 830) deformation space are investigated. It is found that besides the octupole deformation, the triaxiality 
also plays an important role upon the second fission barriers. The non-axial reflection-asymmetric 832 shape in 
some transfermium nuclei with N = 150, namely 46Cm, "Cf 2°°Fm, and ""?No are studied. 


1 Introduction 


“Shape” gives an intuitive understanding of spatial den- 
sity distributions of quantum many-body systems includ- 
ing atomic nuclei. For the description of the nuclear shape, 
it is convenient to use the following parametrization 


4r 


Bau = SARI 


(On), () 
where O,,, are the mass multipole operators. A schematic 
show of some typical shapes is given in Fig. 1 [1]. The 
majority of observed nuclear shapes is of spheroidal form 
which is usually described by 629. Higher-order deforma- 
tions with A > 2 such as £39 also appear in certain atomic 
mass regions [2]. In addition, non-axial shapes in atomic 
nuclei, in particular, the nonaxial-quadrupole (triaxial) de- 
formation >> have been studied both experimentally and 
theoretically [3,4,5]. There is no a priori reason to neglect 
the nonaxial-octupole deformations, especially the 632 de- 
formation [6,7,8]. 

Furthermore, more shape degrees of freedom play im- 
portant roles in the study of potential energy surfaces of 
atomic nuclei. Particularly, various shape degrees of free- 
dom play important and different roles in the occurrence 
and in determining the heights of the inner and outer bar- 
riers in actinide nuclei (in these nuclei double-humped fis- 
sion barriers usually appear). For example, the inner fission 
barrier is usually lowered when the triaxial deformation is 
allowed, while for the outer barrier the reflection asym- 
metric (RA) shape is favored [9, 10, 11, 12, 13]. Nowadays, 
it becomes more and more desirable to have accurate pre- 
dictions of fission barriers also for superheavy nuclei [14, 
15,16, 17,18, 19]. It is usually customary to consider only 
the triaxial and reflection symmetric (RS) shapes for the 
inner barrier and axially symmetric and RA shapes for the 
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Fig. 1. (Color online) A schematic show of some typical nuclear 
shapes: (a) sphere; (b) prolate spheroid; (c) oblate spheroid; (d) 
hexadecapole shape; (e) triaxial ellipsoid; (f) reflection symmet- 
ric octupole shape; (g) tetrahedron; (h) reflection asymmetric oc- 
tupole shape with very large quadrupole deformation. Taken from 
Ref. [1]. 


outer one [15,20,21,22]. The non-axial octupole deforma- 
tions are considered in both the macroscopic-microscopic 
(MM) models [23] and the non-relativistic Hartree-Fock 
theories [24]. 


In order to give a microscopic and self-consistent study 
of the potential energy surface with more shape degrees of 
freedom included, multi-dimensional constrained covari- 
ant density functional theories are developed recently [25, 
26]. In these theories, all shape degrees of freedom 6,, de- 
formations with even yu are allowed, e.g., B20, 822, B30, B32, 
P40, B42, B44, and so on. In this contribution, we present two 
recent applications of these theories: the potential energy 
surfaces of actinide nuclei and the non-axial reflection- 
asymmetric 832 shape in some transfermium nuclei. In Sec- 
tion 2, the formalism of our multi-dimensional constrained 
covariant density functional theories will be given briefly. 
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The results and discussions are presented in Section 3. Fi- 
nally we give a summary in Section 4. 


2 Formalism 


The details of the formalism for covariant density func- 
tional theories can be found in Refs. [27, 28, 29, 30, 31,32]. 
The CDFT functional in our multi-dimensional constrained 
calculations can be one of the following four forms: the 
meson exchange or point-coupling nucleon interactions com- 
bined with the non-linear or density-dependent couplings [25, 
26]. Here we show briefly the one corresponding to the 
non-linear point coupling (NL-PC) interactions. The Start- 
ing point of the relativistic NL-PC density functional is the 
following Lagrangian: 


L£ = V (iy ð" dai Ms) y lai Liin — Ln = Laer = Leou, (2) 
where 
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are the linear, non-linear, and derivative couplings and the 
Coulomb part, respectively. Mg is the nucleon mass, as, 
Av, ATS, ATV, Bs, Vs. VV. Os, Ov, OTS; and OTV are cou- 
pling constants for different channels and e is the electric 
charge. ps, prs, jv, and jry are the iso-scalar density, iso- 
vector density, iso-scalar current, and iso-vector current, 
respectively. The densities and currents are defined as: 


ps = W, (4) 
prs = Vrv, (5) 
dy = vrv, (6) 
Sry = Vrv. (7) 


Starting from the above Lagrangian, using the Slater de- 
terminants as trial wave functions and neglecting the Fock 
term as well as the contributions to the densities and cur- 
rents from the negative energy levels, one can derive the 
equations of motion for the nucleons. Furthermore, for sys- 
tems with time reversal symmetry, only the time-like com- 
ponents of the vector currents (6) and (7) survive. The re- 
sulted equation for the nucleons reads 


hy; =(a- p+ PM 4 S(r)) £ V(r) Vi = evi, (8) 


where the potentials V(r) and S (r) are calculated as 


S = asps + Bsps + YsPS + Os dps 
+(arsprs + 6rsAprs)T3, (9) 
V < avpPv + Yvp$ + ôyApy W 


+ (arvprv + 6rvApry) T3. (10) 


An axially deformed harmonic oscillator (ADHO) ba- 
sis is adopted for solving the Dirac equation [25, 26,33]. 
The ADHO basis are defined as the eigen solutions of the 
Schrodinger equation with an ADHO potential [34,35], 


i 
oa + Vg(z,p)| Palro) = EgPa(ro), (11) 


where 


1 
Vez, P) = 5M(wpp? + w22"), (12) 
is the axially deformed HO potential and w, and w, are the 
oscillator frequencies along and perpendicular to z axis, 
respectively. The solution of Eq. (11) reads 
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Palro) = Can, Ry Pa Ya (0), (13) 
T 
where ġ„,(z) and Ra P) are the HO wave functions, 
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Xs, is a two component spinor and Ce is a complex num- 
ber inserted for convenience. Oscillator lengths b, and b, 
are related to the frequencies by b; = 1/ VMo, and b, = 
1/ (May. 

These basis are also eigen functions of the z component 
of the angular momentum j; with eigen values K = mj+ms. 
For any basis state ®,(ro), the time reversal state is de- 
fined as (ro) = T Palro), where T = ic,,K is the time 
reversal operator and K is the complex conjugation. Appar- 
ently we have Ka = —K, and ma = ^e. These basis form 
a complete set for expanding any two-component spinors. 
For a Dirac spinor with four components, 


_ (Loh? bro) 
Gi (5 fovea) 
where the sum runs over all the possible combination of the 
quantum numbers x = (n,,n,,m;,m,), and f? and gf are 
the expansion coefficients. In practical calculations, one 
should truncate the basis in an effective way [25,26, 33]. 

We expand the potentials V(r) and S(r) and various 
densities in terms of the Fourier series, 


l 
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The nucleus is assumed to be symmetric under the V4 group, 
that is, for all the potentials and densities we have 


Kle, z + 9, 2) = fp, 27 - 9,2) = fp, pz). 


Thus the components f,'s satisfy fa = fa = fa and all the 
terms with odd u vanish. The expansion Eq. (16) can be 
simplified as 


fp,9.D< >) Hoz 


H=- 


exp(iuy). (16) 
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are real functions of p and z. 
The total energy of a nucleus reads 


Etal = [er (x up] (œ - p + BM) We 
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where the center of mass correction E, , can be calculated 
either phenomenologically or microscopically. 

The intrinsic multipole moments are calculated from 
the vector densities by 


O), = f Proy Y uO), (21) 


where Ya, (Q) is the spherical harmonics and 7 refers to the 
proton, neutron or the whole nucleus. 

The potential energy surface (PES) is obtained by the 
constrained self-consistent calculation, 


Ne 
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where H is the RMF Hamiltonian, O,'s are the multipole 
operators to be constrained and N; is the dimension of the 
constraining space. 

Both the BCS approach and the Bogoliubov transfor- 
mation are implemented in our model to take into account 
the pairing effects. For the pairing force, we can use a delta 
force or a separable finite-range pairing force [36, 37,38]. 
More details of the multi-dimensional constraint covariant 
density functional theories can be found in Refs. [25,26]. 


3 Results and discussions 


3.1 One-, two-, and three-dimensional potential 
energy surface of "Pu 


In Ref. [25], one- (1-d), two- (2-d), and three-dimensional 
(3-d) constrained calculations were performed for the ac- 


tinide nucleus Pu. The parameter set PC-PKI is used [39]. 


In Fig. 2 we show the 1-d potential energy curves (PEC) 
from an oblate shape with 8,0 about —0.2 to the fission 
configuration with 829 beyond 2.0 which are obtained from 


E (MeV) 


Fig. 2. (Color online) Potential energy curves of 74°Pu with var- 
ious self-consistent symmetries imposed. The solid black curve 
represents the calculated fission path with V, symmetry imposed: 
the red dashed curve that with axial symmetry (AS) imposed, the 
green dotted curve that with reflection symmetry (RS) imposed, 
the violet dot-dashed line that with both symmetries (AS & RS) 
imposed. The empirical inner (outer) barrier height Bemp is de- 
noted by the grey square (circle). The energy is normalized with 
respect to the binding energy of the ground state. The parameter 
set used is PC-PK1. Taken from Ref. [25]. 


calculations with different self-consistent symmetries im- 
posed: the axial (AS) or triaxial (TS) symmetries combined 
with reflection symmetric (RS) or asymmetric cases. The 
importance of the triaxial deformation on the inner barrier 
and that of the octupole deformation on the outer barrier 
are clearly seen: The triaxial deformation reduces the inner 
barrier height by more than 2 MeV and results in a better 
agreement with the empirical datum; the RA shape is fa- 
vored beyond the fission isomer and lowers very much the 
outer fission barrier. Besides these features, it was found 
for the first time that the outer barrier is also considerably 
lowered by about 1 MeV when the triaxial deformation is 
allowed. In addition, a better reproduction of the empirical 
barrier height can be seen for the outer barrier. It has been 
stressed that this feature can only be found when the axial 
and reflection symmetries are simultaneously broken [25]. 

In order to see how the PES of *4°Pu becomes unstable 
against the triaxial distortion, 2-d PES’s from calculations 
without and with the triaxial deformation were compared 
in Fig. 3 [25]. When the triaxial deformation is allowed, 
the binding energy of **°Pu assumes its lowest possible 
value at each (G29, 830) point. At some points, especially 
those around the two saddle points, non-axial solutions are 
favored than the axial ones. The inner barrier height is low- 
ered by about 2 MeV. About | MeV is gained for the bind- 
ing energy at the second saddle point due to the triaxiality. 
In the regions around the ground state and in the fission 
isomer valleys, only axially symmetric solutions are ob- 
tained. 


A full 3-d PES has been obtained for 74°Pu [25]. In 
Fig. 4 are shown only five typical sections of the 3-d PES of 
40Py in the (822,830) plane calculated around the ground 
state, the first saddle point, the fission isomer, the second 
saddle point, at a point beyond the outer barrier, respec- 
tively. The following conclusions were drawn by examin- 
ing these 3-d PES’s [25]: (1) The ground state and the fis- 
sion isomer are both axially and reflection symmetric. The 
stiffness of the fission isomer is much larger than that of 
the ground state against both the 622 and 630 distortions. 


201707.00938v1 


chinaXiv 


EPJ Web of Conferences 


w/o y 


240 
Pu 


| PC-PK1 


30 


Boo 


Fig. 3. (Color online) Potential energy surfaces of *“°Pu in the 
(820, 830) plane from calculations (a) without and (b) with the tri- 
axial deformation included. The energy is normalized with re- 
spect to the binding energy of the ground state. The numbers in 
(b) show the values of 82 at these points. The fission path is rep- 
resented by a dash-dotted line. The ground state and fission iso- 
mer are denoted by full and open circles. The first and second 
saddle points are denoted by full and open triangles. The contour 
interval is 1 MeV. Taken from Ref. [25]. 
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Fig. 4. (Color online) Sections of the three-dimensional PES of 
40Pu in the (827,839) plane calculated at B29= 0.3 (around the 
ground state), 0.6 (around the first saddle point), 0.9 (around the 
fission isomer), 1.3 (around the second saddle point) and 1.6 (be- 
yond the outer barrier), respectively. The energy is normalized 
with respect to the binding energy of the ground state. The con- 
tour interval is 0.5 MeV. Local minima are denoted by crosses. 
Taken from Ref. [25]. 
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Fig. 5. (Color online) The inner (Bi) and outer (Bf) barrier heights 
of even-even actinide nuclei. The axial (triaxial) results are de- 
noted by open (full) symbols. The empirical values are taken from 
Ref. [40] and represented by grey sguares. Taken from Ref. [25]. 


(2) The second saddle point which is close to B29 = 1.3 ap- 
pears as both triaxial and reflection asymmetric shape. (3) 
The triaxial distortion appears only on the top of the fission 
barriers. 

From the investigation of the one-, two-, and three- 
dimensional PES of "Pu, we can learn a lot about the 
importance of different shape degrees of freedom in differ- 
ent regions of PES in actinide nuclei. These information 
could be useful in further systematic calculations. 


3.2 Inner and outer fission barriers of even-even 
actinide nuclei 


Guided by the features found in the 1-d, 2-d, and 3-d PES’s 
of Pu, the fission barrier heights were extracted for even- 
even actinide nuclei [25]. The calculated values were com- 
pared with empirical ones recommended in RIPL-3 [40]. 

As it has been shown previously, around the inner bar- 
rier an actinide nucleus assumes triaxial and reflection sym- 
metric shapes. Thus in order to obtain the inner fission 
barrier height we can safely make a one-dimensional con- 
straint calculation with the triaxial deformation allowed 
and the reflection symmetry imposed. In Fig. 5(a) we show 
the calculated inner barrier heights Bi and compare them 
with the empirical values. It is seen that the triaxiality low- 
ers the inner barrier heights of these actinide nuclei by 
1 ~ 4 MeV as what has been shown in Ref. [13]. In general 
the agreement of our calculation results with the empirical 
ones is very good with exceptions in the two thorium iso- 
topes and **°U. Possible reasons for these disagreements 
were discussed in Ref. [25]. 
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To obtain the outer fission barrier height Bẹ, the situ- 


ation becomes more complicated because more shape de- 
grees of freedom play important roles around the outer fis- 
sion barrier. In Ref. [25], 2-d constraint calculations were 
made carefully around the second saddle points for even- 
even actinide nuclei. In the lower panel of Fig. 5 we show 
the results of outer barrier heights B? and compare them 
with empirical values. For most of the nuclei investigated 
here, the triaxiality lowers the outer barrier by 0.5 ~ 1 
MeV, accounting for about 10 ~ 20% of the barrier height. 
One finds that the calculation with the triaxiality agrees 
well with the empirical values and the only exception is 
248Cm. From the calculation with the axial symmetry im- 
posed, the outer barrier height of "Cm is already smaller 
than the empirical value. The reason for this discrepancy 
may be related to that there are two possible fission paths 
beyond the first barrier [25]. 

In Ref. [25], it was also examined the parameter de- 
pendency of the influence of triaxiality on the outer fission 
barrier and the lowering effect of the triaxiality on the outer 
fission barrier was also observed when parameter sets other 
than PC-PKI are used. 


3.3 Non-axial octupole shapes in N = 150 isotones 


Nowadays the study of nuclei with Z ~ 100 becomes more 
and more important because it not only reveals the struc- 
ture of these nuclei themselves but also provides significant 
information for superheavy nuclei [41,42,43,44]. One of 
the relevant and interesting topics is how to explain the 
low-lying 2~ states in some N = 150 even-even nuclei. In 
these nuclei, the bandhead energy E(2~) of the lowest 27 
bands is very low [45]. It is well accepted that the octupole 
correlation is responsible for it. For example, a quasiparti- 
cle phonon model with octupole correlations included was 
used to explain the excitation energy of the 2" state of 
the isotones with N = 150 [46]. In Ref. [47], Chen et 
al. proposed that the non-axial octupole Y32-correlation re- 
sults in the experimentally observed low-energy 27 bands 
in the N = 150 isotones and the reflection asymmetric shell 
model calculations reproduces well the experimental ob- 
servables of these 2~ bands. 

The non-axial reflection-asymmetric 832 shape in some 
transfermium nuclei with N = 150, namely 246Cm, %8Cf, 
250Fm, and ??No were investigated with the multidimen- 
sional constrained covariant density functional theory [48]. 
The parameter set DD-PC1 is used [49]. For the ground 
states of “Cf and °°Fm, the non-axial octupole deforma- 
tion parameter 632 > 0.03 and the energy gain due to the 
32 distortion is larger than 300 keV. In *4°Cm and **7No, 
shallow 632 minima are found. 

The triaxial octupole Y32 effects stem from the cou- 
pling between pairs of single-particle orbits with A j = Al = 
3 and AK = 2 where j and / are the total and orbit angular 
momenta of single particles respectively and K the projec- 
tion of j on the z axis. In Fig. 6, we show the proton and 
neutron single-particle levels near the Fermi surface for 
48 Cf as a function of 32 with B20 fixed at 0.3. It was shown 
that the spherical proton orbitals 72 f7/2 and 71i}3/2 are very 
close to each other [48]. This near degeneracy results in oc- 
tupole correlations. As seen in Fig. 6, the two proton levels, 
[521]3/2 originating from 27/2 and [633]7/2 originating 
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Fig. 6. (Color online) The single-particle levels near the Fermi 
surface for (a) protons and (b) neutrons of "Cf as a function of 
B32 with B20 fixed at 0.3. 


from 1i)3/2, satisfying the 4j = Al = 3 and AK = 2 con- 
dition, are very close to each other at 629 = 0.3. Therefore 
the non-axial octupole Y32 develops and with 832 increas- 
ing from zero, an energy gap appears at Z = 98. Similarly, 
the spherical neutron orbitals v2g9,> and vl j;s,> are very 
close to each other [48]. The neutron levels [734]9/2 origi- 
nating from 1 j15;2 and [622]5/2 originating from 2g9;2 are 
also close to each other and they just lie above and be- 
low the Fermi surface. This leads to the development of a 
gap at N = 150 with 832 increasing. The Y32 correlation in 
N = 150 isotones is from both the proton and the neutron 
and for “8Cf the correlation is the most pronounced [48]. 


4 Summary 


In this contribution we present some applications of the 
multi-dimensional constrained covariant density functional 


EPJ Web of Conferences 


theories in which all shape degrees of freedom 6,, defor- 
mations with even yu are allowed. The potential energy sur- 
faces of actinide nuclei in the (G29, 822, 830) deformation 
space are investigated. It is found that besides the octupole 
deformation, the triaxiality also plays an important role 
upon the second fission barriers. For most of even-even 
actinide nuclei, the triaxiality lowers the outer barrier by 
0.5 ~ 1 MeV, accounting for about 10 ~ 20% of the barrier 
height. The non-axial reflection-asymmetric 832 shape in 
some transfermium nuclei with N = 150, namely 246Cm, 
248Cf, 250Em, and ??No are studied. Due to the interac- 
tion between a pair of neutron orbitals, [734]9/2 originat- 
ing from vj15/2 and [622]5/2 originating from vgo/2, and 
that of a pair of proton orbitals, [521]3/2 originating from 
zfin and [633]7/2 originating from 7i}3/2, rather strong 
non-axial octupole Y32 effects have been found for 248 Cf 
and *°Fm which are both well deformed with large axial- 
quadrupole deformations, 629 ~ 0.3. For 246Cm and "No, 
a shallow minima develops along the 632 deformation de- 
gree of freedom. 
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